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Abstract 



The effect of ring vortices on the destabilization of supersaturated liquid 
helium mixtures at very low temperatures is investigated as a function of pres- 
sure. We have found that large ring vortices trigger the segregation of 3 He at 
smaller values of the concentration than small ring vortices. Our calculations 
indicate that the existence of ring vortices in the mixture is a possible mech- 
anism to understand the rather small degree of critical supersaturation found 
in the experiments. 
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I. INTRODUCTION 



The study of supersaturated liquid helium mixtures at very low temperatures (T) has 
received a renewed experimental interest [0-|v|. The main goals of these experiments have 
been, on the one hand, to determine the crossover from the classical thermal regime, to the 
quantum regime in the nucleation process of the pure 3 He phase from the 3 He- 4 He mixture, 
and on the other hand, to study the mixture in the metastable supersaturated state, whose 
existence has been established long ago ||0]. 

One of the more interesting observations of recent experiments is the rather small degree 
of supersaturation achieved in them, of the order of 1 % in [0,|4]|| (see also ||), and below 
0.5 % in . Similar supersaturations were found in [00], but at that time no systematic 

studies of the metastable phase were undertaken. 

An extrapolation of the measured 3 He chemical potential excess Afi% along the coex- 
istence line yielded dA/i^/dx > up to values of the 3 He concentration x > 16 % [ 1Q| , 
indicating that the mixture could be in the metastable, homogeneous phase up to that con- 
centration. That belief seemed further sustained by the calculations of Lifshitz et al in their 



classical paper about the kinetics of the nucleation process in liquid helium mixtures |LT 



and by microscopic JT2| and density functional [HJ estimates of the spinodal line in the P — x 
plane at zero temperature. 

Clearly, there is nothing a priori contradictory between the results of fI"2| , [T3|j and the 
small critical supersaturation degree Ax cr found experimentally. Likely, the segregation of 
3 He from the mixture occurs before reaching the spinodal line. What is really intriguing is 
that nucleation calculations [P, p] , p^1 , which incorporate the basic physics to consider 3 He 
droplets as nucleation seeds for phase separation, have yielded results that are around one 
order of magnitude larger than experiment. Other calculations have put the emphasis in the 
diffusion mechanism of the quantum decay of the metastable mixture and in the quantum- 



to-thermal crossover temperature T* |fL5|| , taking as an input the experimental Ax cr . Yet, 
T* comes out a factor of four larger than experiment. 

It is then quite natural to look for another physical mechanism as a possible explanation 
for the disagreement. One such mechanism is heterogeneous nucleation. However, the system 
is extremely pure, and the experimental cells are covered with 4 He, avoiding nucleation on the 
walls. Since at the low temperatures involved in these experiments, the solubility of 3 He into 
4 He is limited and 4 He is still superfluid, one is thus left with an obvious candidate, namely 
heterogeneous nucleation on vortices whose core is filled with 3 He [ 16| . This possibility was 
already mentioned in ||, independently put forward and worked out in same detail in |14 
and further developed in |17]], as well as in |fL8| . 



The basis of the model employed in |14| is that a 4 He vortex line described by a hollow 
core model (HCM) is gradually filled with 3 He as x increases. In the metastable region, 
these lines become unstable above a certain Ax value, and the resulting Ax cr turns out 
to be in the right order of magnitude (around 1.5 %), especially as compared with the 
experimental results of [0|J. The HCM results have been confirmed by a more elaborated 
density functional calculation |T7 . 



The purpose of the present work is to complete the above scenario, studying the effect 
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that a more realistic vortex geometry, namely a ring vortex, may have on the quantitative 
results, and also to compare them with the recent experimental results of ||. As in previous 
works [p^|TT|,|T8l , the presence of vortices is taken for granted. Indeed, it has been recognized 



for a long time that 'any container of superffuid 4 He, treated in conventional fashion, will be 
permeated ab initio by numerous quantized vortices stabilized by surface pinning' |19| . The 
interesting problem of vorticity nucleation in liquid 4 He, and the effect of 3 He impurities on 
it has been addressed in detail in the past (see for example [ffi|-|23|). 

This work is organized as follows. We present in Section II the model we have used to 
describe the evolution of a ring vortex as a function of pressure and 3 He concentration. The 
numerical results are collected in Section III. Although our main interest is in the metastable 
region in the pressure-concentration (P — x) plane at very low temperatures, we have also 
paid some attention to the stable region, in view that previous studies have not dealt with 
appreciable 3 He concentrations. A discussion is presented in Section IV, and some technical 
details are given in an Appendix. 



II. RING VORTEX ENERGETICS 

To carry out the program we have outlined in the introduction, we first determine the 
kinetic energy of the ring vortex (RV), for which we assume a hollow core model, i.e., the 
4 He density is zero within the vortex core, and equal to the bulk value at the given pressure 
elsewhere. Following [^,^(J], we consider the liquid to be incompressible (V • v ~ 0), so 
that one can introduce a velocity vector potential A(r ) such that v = V x A. The vector 
potential is fixed if the vorticity uj = V x v is specified: 

w = Vx(Vxi) . (1) 

We make use of the line-source approximation for the vorticity. For a ring vortex with 
radius of curvature R, circulation number n = 1, 2, and quantum circulation ko = nh/m^ 
located on the z = plane, it means that 

u(r) = k 5(r-R)5(z)4> , (2) 

where (r, z, 0) are the cylindrical coordinates and <fi = (— sin(f), coscj), 0) is a unit vector. 
Unless explicitly indicated, we shall be considering n = 1 as circulation number. It is 
straightforward to show that the integral solution of Eq. (U) reads A(f) = A (r, z)4>, with 

A (r z) = —R cos(p'd(p' 

0[r,Z> An Jo ^/r' 1 + R 1 -2rRcos(j)' + z 2 ' U 

The velocity of the superffuid outside the RV can be obtained as 

dA 

n '- = — 

az 

(4) 

r or 
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and the kinetic energy of the superfluid is evaluated as an integral over the core surface S 
of the ring vortex: 



777*4 f 

E k in = — —Pa / A (r, z) (v z f - v r z) ■ hdS 
l Js 



(5) 



where the normal n is directed into the 4 He fluid. The boundary condition that no velocity 
line crosses the vortex core imposes that v ■ n — and yields 



n 



| v | | v | 

as well as the differential equation obeyed by the velocity lines: 

dr v r 

dz v z 



(6) 



(?) 



Using Eq. (]4j), it is easy to verify that solving Eq. ([/]) to obtain the velocity lines is 
equivalent to solving the algebraic equation 



rA (r, z) = Constant . 



(8) 



We have specified the cross section profile of the vortex core, compatible with the shape of 
the velocity lines, introducing a length a such that (R — a) is the closest distance of the core 
surface to the origin of coordinates: 



(R-a) A (R-a,0) = C a 



(9) 



That fixes the constant in Eq. (|) and permits one to obtain the core profile. In the case 
of large ring vortices with R>> a, for which a circular torus is a good approximation, the 
length a represents the radius of the circular cross section of the core. We have numerically 
checked that if a/R « 1, A at the surface of the core, and the kinetic energy Ekin are 
indeed those of the hollow RV given for example in 



A = ^ 
2vr 



E, 



m 4 



Lin 



' 8R 
In — 

a 



(10) 



In 



8R 



We give in the Appendix the general expressions of A , v r and v z which are needed to 
evaluate the kinetic energy, Eq. (H). 

Figure |T] shows several cross sections of the core for R = 30 A and a = 10, 15 and 20 
A, respectively. One can see that these cross sections are very non-circular, especially when 
the ratio a/R is large. 

It is now a simple task to write the energy gain caused by the presence of a ring vortex in 
the homogeneous mixture at given P. As x increases, the 3 He atoms located at the surface 
of the vortex migrate to the interior of the hollow core. Within the HCM, one has [O] 



E v (R,a) = aS-A f i 3 p 3 V + E i 



kin 



:n) 
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where V is the volume of the core, p 3 the 3 He atom density inside the vortex core, a the 
surface tension of the mixture-pure 3 He interface, and A/i 3 is the 3 He chemical potential 
excess, i.e., the difference between the chemical potential of 3 He in the metastable mixture, 
and in pure 3 He. Consequently, A/i 3 has a sign. It is negative when the mixture is stable, 
and positive when it is not. 

The reliability of eq. (|TT|) to describe 4 He vortices filled with 3 He has been checked in 



T7[ for the case of vortex lines using a density functional approach which takes into account 
the surface diffuseness of both liquids at the interface, and the kinetic energy of the 4 He 
atoms that penetrate into the vortex core. Altogether, both effects have been found to be 
small, the reason being that the core radii of critical vortices are large, of the order of 10 
A, and consequently, the sharp surface approximation embodied in the HCM works well. 
For example, at zero T and P, the HCM yields Ax cr ~ 1.6 % [14]], whereas the density 
functional theory yields Ax cr ~ 1.8 % |L7|. 



The temperatures of interest here are below 150 mK ( Bearing in mind that 



we are looking for a sizeable decrease in the value of Ax cr with respect to that of |]TT 
and that thermal effects are not expected to play an important role because of the small 
T-dependence of o ]27| , A/i 3 |L0| and p 3 p8[ in that temperature domain, we shall carry out 
a T = calculation. 

At given P and 3 He concentration x = p 3 / (p 3 + p4), where p 4 is the 4 He atom density, the 
thermodynamical characteristics of the liquid have been described by a density functional 
that among other properties, reproduces accurately the maximum solubility line and surface 
tension of the interface as a function of P ||29|| . We give in Fig. ^| Ap 3 as a function of x 
for P = 0, 5 and 10 atm. The Ap 3 curves intersect the x-axis at a concentration x s (P) 
that corresponds to the maximum solubility of 3 He into 4 He for the given P. The slopes 
<9(Ap 3 )/<9x \ Xs will be used in the next section to determine the critical supersaturation of 



the mixture when ring vortices are present, transforming the Ap 3 values used in Eq. (|TT|) 
into Ax values as 



Ax ~ Ap 3 



dx 



(12) 



This is justified in view of the smallness of Ax. We have obtained <9Ap 3 /<9x \ Xg = 2.49, 1.97 
and 1.85 K, and x s = 6.40, 9.01 and 9.58 % for P— 0, 5 and 10 atm, respectively. 



III. RESULTS 

A. Stable phase-diagram region 

We have first considered the region in the P — x plane where the homogeneous mixture 
is stable, i.e., P > and Ax = x — x s < 0, and have minimized the energy Eq. ( pT|) of a 
RV of radius of curvature R with respect to a. The value of R is considered as an intrinsic 
characteristic of the vortex that depends on the process of formation, whereas the a value 
may change with changing x, and its optimum value for given x has to be found through 
minimization. 



As expected from our previous works on vortex lines [14,17 , in this case there is only one 



extremum of E V (R, a) corresponding to a stable ring vortex. The energy of the equilibrium 
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configuration and the effective radius of the vortex core do, defined as nal = S, are drawn 
in Figs. H and [| as functions of R for P = and 10 atm, and Ax = -1, -3 and -5 %. 



B. Metastable phase-diagram region 

When Ax is positive, E v (R,a) at fixed R has in general two extrema, one maximum 
and one minimum: the ring vortex is no longer stable and may trigger the demixing process 
of 3 He if it overcomes the barrier whose height AE V is the energy difference between both 
extrema. It may happen either by thermal or by quantum fluctuations. One such energy 
barrier is drawn in Fig. |5] for an R = 100 A ring vortex at P = and Ax = 1 %. We show 
in Fig. P the barrier heights as a function of Ax, for different R and P values. 

The barrier height decreases with increasing Ax, and there is a critical value Ax cr at 
which it disappears. This defines a critical ring vortex that depends on the values of R and 
P, and may freely expand since it is barrierless. 

We display in Fig. |7| the effective radius of the critical vortex as a function of R for 
P = 0, 5 and 10 atm. We see that in contradistinction with the stable case, the core radii 
are large, and as in the vortex line case, we expect the sharp surface approximation to work 
well. The critical configurations corresponding to R = 30 A and P = 0, 5 and 10 atm are 
shown in Fig. §. The associated Ax cr values are 2.75, 1.82 and 1.37 %, respectively. For R 
= 100 A, these values are 1.67, 1.04 and 0.76 %. 

From the previous discussion, it is clear that one can define, for a given i?-vortex, a kind 
of ring vortex spinodal line in the P — x plane, lying in the metastable region, beyond which 
a liquid helium mixture hosting i?-vortices is absolutely unstable. Such spinodal lines can 
be drawn combining the appropiated (P, Ax) values from Fig. |6|. We have represented in 
Fig ^| Ax cr as a function of R for P = 0, 5 and 10 atm. From this figure one can see that 
fixed P, large ring vortices have associated smaller values of Ax cr . It means that large ring 
vortices trigger the seggregation of 3 He at smaller supersaturations that small ring vortices. 



IV. DISCUSSION 

Very recently, Chavogets et al || have carried out a systematic investigation of nucleation 
in superfluid 3 He- 4 He solutions at T = 125 mK as a function of P, which is especially well 
suited for a detailed comparison with the present calculations since the temperature is kept 
constant at a value of 125 ± 5 mK. However, these authors stress that it cannot be excluded 
that the experimental results are influenced by the way the phase transition is initiated, 
since different methods are used to attain the initial metastable state. 



We show in Fig. 10 the experimental results of || for Ax cr (dots) and the spinodal 
values (solid lines) corresponding to ring vortices of R — 100, 200 and 500 A. For R = 30 
A the Ax cr values go from ~ 2.7 % at P = to ~ 1.5 % at P = 10 atm. The calculated 
Ax cr 's tend to decrease with P, but still are a factor of two larger than the experimental 
ones even for the most favorable case of large ring vortices. However, the improvement over 
nucleation calculations that do not consider the possible existence of vortices in the mixture 
is apparent P,|8|,[TT| . 
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The major discrepancy appears at pressures close to zero atm. Chagovets et al || have 
found that Ax cr has a maximum of 1.5 % around P = 0.1 atm, and that Ax cr drops to ~ 
0.7 % around P = 0.02 atm. We have not found this behavior. However, we would like 
to recall that at pressures close to zero, the emergence of a free liquid-vapor interface may 
help nucleate the new phase from the 3 He-enriched surface layer [Q . That might explain the 
small Ax cr found at low pressures, although the authors of || have discarded this possibility 



The agreement between theory and experiment can be dramatically improved considering 
the existence of ring vortices with circulation number n > 1. These vortices are rare, since 
creating an i?-vortex with n = 2 costs more energy than creating two .R-vortices with n = 
1. This is the reason why in classical experiments where vortices are nucleated and dragged 
by moving ions only n = 1 vortices are detected Still, it is interesting to look at their 



effect on Ax cr . Roughly, one should expect that large RV with n = 2 originate Ax cr four 
times smaller than vortices with n = 1 (see for example the discussion after Eq. (4) of [0), 
thus bringing theory closer to experiment. This is indeed the it can be seen from the 

dashed lines in Fig. [10] which represent the spinodal values corresponding to ring vortices 
with n = 2 and R = 30, 50, 75, and 100 A. 

The above comparisons have involved the calculated spinodal Ax cr values, which are 
an upper bound to those one would obtain considering thermal activation over the energy 
barrier (quantum tunneling through it is only effective for temperatures below a few tens 
of niK, see We have checked that, although possible, thermal activation plays a minor 
role. Let us consider, as an example, the case of the R = 50, n = 2 ring vortex at T = 125 
niK, P = 10 atm, for which the spinodal value is Ax CT = 0.56 %. Within classical nucleation 
theory, the decay rate T, i.e., the number of RV 'passing' over the barrier per unit time and 
volume due to thermal fluctuations is given by 

r = T exp{-AE v /kT) , (13) 

where the prefactor T can be roughly estimated as 

r ° = m ■ (14) 

Solving the equation 1 = V exp t exp T, taking as experimental volume V exp and time t exp 
those of ||, and relating the AE V value to Ax, one gets Ax cr = 0.54 % for AE V = 10 K, 
and Ax cr = 0.55 % for AE V = 8 K, which are very close to the spinodal values. A similar 



conclusion was drawn for vortex lines in 14 



Quantum decay is effective only below a crossover temperature T* ~ 10 mK ||. None 
of the calculations carried out so far is able to reproduce Ax cr and T* simultaneously, and 
the problem is still open. In the case of vortex lines, T* has been calculated to be ~ 0.1 mK 



in ||18|| . Using a functional integral method similar to that of ||32j| , we have obtained T* ~ 
1 mK. The smallness of this value stems from the fact that in the calculations, the decay 
process is always found to take place very close to the spinodal line. 
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APPENDIX: 



In this Appendix we give the expressions and method we have used to numerically obtain 
the velocity lines. From Eqs. (|3|4J) one has: 

r*/a 2 cos 2 6 - 1 



v r (r,z) 



k Rz 



(Al) 



vr [(r + R) 2 + z 2 f 2 Jo (1 - s cos 2 ^/ 2 

where we have defined s = 4rR/ [(r + R) 2 + z 2 ]. The integral Eq. ( |A1| ) is written in terms 
of hypergeometric functions |33j yielding: 

fco Rz 



v r (r,z) 
Similarly, 

v z (r, z) 



2 [(r + R) 2 + z 2 f 2 
k Q R 



[F(3/2,3/2;2;s)-F(l/2,3/2;l;s)] 



(A2) 



+ 



2r[(r + R) 2 + z 2 ^ /2 
k R 



m [F(3/2, 1/2; 2; s) - F(l/2, 1/2; 1; s)\ 



2 [(r + i?) 2 + z 2 f 2 
3R 

+ — F(5/2,3/2;3;s)] 



[(r + i?) F(l/2, 3/2; 1; s) - (r + 2/2) F(3/2, 3/2; 2; s) 



and 



Aa(r,z) 



k R 



^[F(3/2,l/2;2;s)-F(l/2,l/2;l;s)] . 



(A3) 



(A4) 



2[(r + R) 2 + z 2 } 1/2 

The hypergeometric functions entering Eqs. (|A2HA4j) can be written in terms of the 



complete elliptic integrals E and K [S3J as follows: 



F(l/2,l/2;l;s) = -K(V5) 

7T 

F(l/2, 3/2; 1; s) = F(3/2, 1/2; 1; s) = 1^3. 

7T 1 — S 



F(3/2, 1/2; 2; s) = F(l/2, 3/2; 2; s) = — [K(^) - E(^)] 

7TS 



F(3/2,3/2;2;s) = — 

7TS 



E(v^) 



F(5/2,3/2;3;s) 



16 

37TS 5 



1 - s 
2 - s 



K(v^) 



1 - s 



E(v^) -2K(Vs) 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 



To evaluate E and K we have used accurate polynomial approximations |34 
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FIGURES 



FIG. 1. Core cross sections corresponding to ring vortices of radius of curvature R = 30 A 
and a-values of 10, 15 and 20 A. The vorticity line is indicated by a dot. 

FIG. 2. 3 He chemical potential excess [/^(x) — ^3 {pure)} (K) as a function of x (%) for P = 
0, 5 and 10 atm. 

FIG. 3. Bottom panel: Ring vortex energy E v (K) function of R (A) for P = 0. Top 
panel: same as bottom panel for P = 10 atm. Both panels correspond to stable vortices, and the 
curves represent the results for Ax = -1 % (solid line), -3 % (dotted line) and -5 % (dashed line). 

FIG. 4. Bottom panel: Effective radius ao (A) of the vortex cross section as a function of R 
(A) for P = 0. Top panel: same as bottom panel for P = 10 atm. Both panels correspond to 
stable vortices, and the curves represent the results for Ax = -1 % (solid line), -3 % (dotted line) 
and -5 % (dashed line). 

FIG. 5. Ring vortex energy E v (K) as a function of a (A) for an R = 100 A ring vortex at P 
= and Ax = 1 %. 

FIG. 6. Barrier heights AE V (K) as a function of Ax for P = 0, 5 and 10 atm. From top to 
bottom, the panels correspond to R = 10, 30 and 50 A. 

FIG. 7. Effective radius ao (A) of the critical vortex as a function of R (A) for P = 0, 5 and 
10 atm. 

FIG. 8. Core cross sections of critical ring vortices corresponding to R = 30 A and P = 0, 5 
and 10 atm. 

FIG. 9. Critical supersaturation Ax cr (%) as a function of R for P = 0, 5 and 10 atm. 

FIG. 10. Critical supersaturation Ax cr (%) as a function of P (atm). The dots are the 
experimental values of Ref. [8]. From top to bottom, the solid lines represent the spinodal values 
for ring vortices with R = 100, 200 and 500 A and circulation number n = 1, and the dashed 
lines represent the spinodal values for ring vortices with R = 30, 50, 75 and 100 A and circulation 
number n = 2. 
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